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We investigate the effects of chiral three-nucleon force (3NF) on proton scattering at 65 MeV and ^He scat¬ 
tering at 72 MeV/nucleon from heavier targets, using the standard microscopic framework composed of the 
Brueckner-Hartree-Fock (BFIF) method and the ^-matrix folding model. For nuclear matter, the g matrix is 
evaluated from chiral two-nucleon force (2NF) of N^LO and chiral 3NF of NNLO by using the BHF method. 

Since the g matrix thus obtained is numerical and nonlocal, an optimum local form is determined from the on- 
shell and near-on-shell components of g matrix that are important for elastic scattering. For elastic scattering, 
the optical potentials are calculated by folding the local chiral g matrix with projectile and target densities. This 
microscopic framework reproduces the experimental data without introducing any adjustable parameter. Chiral- 
3NF effects are small for proton scattering, but sizable for "^He scattering at middle angles where the data are 
available. Chiral 3NF, mainly in the 27r-exchange diagram, makes the folding potential less attractive and more 
absorptive for all the scattering. 

PACS numbers: 21.30.Fe, 24.10.Ht, 25.40.Cm, 25.55.Ci 


I. INTRODUCTION 

An important current issue in nuclear physics is to under¬ 
stand the effects of three-nucleon force (3NF) on finite nuclei, 
nuclear reactions and nuclear matter. Quantitatively decisive 
roles of 3NFs have been established in properties of light nu¬ 
clei as well as of nuclear matter [1]. This issue started with 
the 27r-exchange 3NF proposed by Fujita and Miyazawa [2]. 
Recently, a major breakthrough was made on this issue with 
chiral effective field theory (Ch-EFT); see Refs. [3, 4] and 
references therein. Ch-EFT is a theory based on chiral per¬ 
turbation theory to provide a low-momentum expansion of 
two-nucleon force (2NF), 3NF and many-nucleon forces. Us¬ 
ing this theory, one can define multi-nucleon forces system¬ 
atically. The effects of chiral 3NF were analyzed in many 
papers; e.g., see Ref. [5] for light nuclei, Refs. [6, 7] for 
ab initio nuclear-structure calculations in lighter nuclei and 
Refs. [8-12] for nuclear matter. Recently the role of chiral 
four-nucleon forces was also analyzed for nuclear matter [13]. 
When the g matrix (the effective nucleon-nucleon interaction 
in nuclear medium) is calculated from chiral 2NFH-3NF with 
the Brueckner-Hartree-Fock (BHF) method, it well accounts 
for the empirical properties of symmetric nuclear matter [9]. 
The g matrix depends on the nuclear-matter density p; i.e., 
g = g{p). Chiral-3NF effects become more important as p 
increases. 

Another important issue in nuclear physics is microscopic 
understanding of nucleon-nucleus (NA) and nucleus-nucleus 
(AA) optical potentials. The optical potentials are essential in 
describing not only elastic scattering but also inelastic scat¬ 
tering and transfer and breakup reactions. In fact, the optical 
potentials are essential inputs in distorted-wave Born approx¬ 
imation and continuum discretized coupled-channel method 
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(CDCC) calculations [14]. 

The p-matrix folding model is a standard method of calcu¬ 
lating the optical potential microscopically. Actually, many 
nuclear reactions have been analyzed with the model. In the 
model, the potential is obtained by folding g{p) with projec¬ 
tile and target densities (pp and px) for AA scattering and 
with pt for NA scattering; e.g., see Refs. [15-19] for g{p) 
and Refs. [20-22] for the folding procedure. The model is 
called the single-folding (SF) model for NA scattering and the 
double-folding (DF) model for AA scattering. 

For NA elastic scattering, the SF model based on the Mel¬ 
bourne g matrix [16], constructed from the Bonn-B nucleon- 
nucleon (NN) interaction [23], well reproduces the experi¬ 
mental data with no adjustable parameter. In the folding pro¬ 
cedure, the value of p in p(p) is assumed to be a value of px 
at the midpoint r^a of interacting two nucleons: p = px(fm)- 
Target-excitation effects on the elastic scattering are thus well 
described by the SF model based on the local-density approx¬ 
imation. 

We have recently investigated chiral-3NF effects on the 
description of NA scattering [24] and A A scattering [25], 
r2c_|_i2c and by modifying the Melbourne g ma¬ 

trix by introducing spin- and isospin-dependent multiplicative 
factors to simulate the 3NF effects. Results have shown that 
the 3NF effects are small for NA scattering, because the pro¬ 
cess is mainly governed by the interaction at low density re¬ 
gions. This reinforces the success of calculations using the 
Melbourne g matrix without considering 3NFs. On the other 
hand, sizable effects are found for AA scattering at around 
80 MeV/nucleon through the repulsive contribution in the real 
part and the enhanced absorptive potential. However, these 
are rather exploratory, relying on the Melbourne g matrix. In 
this paper, we present full chiral g matrices parameterized in 
a 3-range Gaussian form on the basis of nuclear matter g- 
matrix calculations with the 2NF and 3NF of Ch-EFT, and 
apply them to the SF and DF models. 

The p-matrix DF model for AA scattering has a practical 
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problem. In nuclear matter calculations, we have to consider 
two Fermi spheres and the g matrix should be obtained by 
solving scattering between a nucleon in a Fermi sphere and 
a nucleon in another Fermi sphere [26, 27], but it is quite 
difficult in practice. In fact, the g matrix is evaluated by 
solving nucleon scattering from a single Fermi sphere. For 
consistency with the nuclear-matter calculation, we assumed 
P = PT(fm) in g{p) and applied this framework to scat¬ 
tering from heavier targets in a wide range of incident ener¬ 
gies from 30 MeV/nucleon to 180 MeV/nucleon [22, 28]. The 
Melbourne p-matrix DF model with the target-density approx¬ 
imation (TDA) well reproduces the data with no adjustable 
parameter, particularly for total reaction cross sections ctr and 
forward differential cross sections. The DF-TDA model does 
not include projectile-excitation effects, but it was confirmed 
by CDCC calculations that the effects are negligible for ^He 
scattering. Precisely, the effects are appreciable at incident en¬ 
ergies lower than 40 MeV/nucleon, but they enhance ctr only 
by a few percent. It is very likely that projectile-excitation ef¬ 
fects are even smaller for ^He scattering, since ^He is less 
fragile than ^He. The practical problem is thus solved for 
^’^He scattering. Therefore, this DF-TDA model is used in 
this paper. 

For heavier projectiles than ^He, it is quite difficult to 
include all projectile-excitation effects explicitly. For such 
AA scattering, the frozen-density approximation (FDA), p = 
Pp('fm) + PT('fm), is often taken as a value of p in g{p), 
although g{p) is obtained by solving nucleon scattering 
on a single Fermi sphere. The DF-FDA model includes 
projectile-excitation effects approximately. The model based 
on the Melbourne g matrix well reproduces measured ctr for 
12 ,14-16(3 [29, 30] and Ne and Mg isotopes [21, 31]. As an 
important result, the microscopic analyses conclude that ^^Ne 
and ^^Mg are deformed halo nuclei. For ^’^He scattering, 
however, the DF-TDA model always yields better agreement 
with the experimental data than the DF-FDA model [22, 28]. 

The g matrix obtained is quite inconvenient in many appli¬ 
cations, since it is nonlocal and numerical. The Melbourne 
group showed that elastic scattering are mainly determined by 
the on-shell part of g{p) [16]. Making a fitting to the on- 
shell and near-on-shell components of the g matrix, the group 
provided g{p) with a local (Yukawa) form in order to make 
the folding procedure feasible [16, 32, 33]. The Melbourne g 
matrix thus obtained accounts for NN scattering in the limit 
of p = 0, and the SF model based on the Melbourne g ma¬ 
trix explains NA scattering systematically with no adjustable 
parameter, as mentioned above. 

In this paper, we consider heavier targets such as "‘^Ca, ®®Ni 
and ^°®Pb to make our discussion clear, since the g matrix is 
evaluated in nuclear matter and the p-matrix folding model 
is considered to be more reliable for heavier targets. Taking 
the Melbourne-group procedure [16, 32, 33], we provide the 
chiral g matrix with a 3-range Gaussian form for each of the 
central, spin-orbit and tensor components, since the Gaussian 
form is much more convenient than the Yukawa form in many 
applications whereas the two forms yield the same results for 
NA and AA scattering. The ranges and the depths of individ¬ 
ual components are determined for each energy and density 


so as to reproduce the on-shell and near-on-shell matrix ele¬ 
ments of the original g matrix. For the central part of g matrix, 
the present ranges of 3-range Gaussian form are (0.4,0.9, 2.5) 
fm and close to those of Ref. [17]. We call the analytic form 
“Gaussian chiral g matrix” and the original numerical g ma¬ 
trix “original chiral g matrix”, when we need to identify the 
two. The folding model based on Gaussian chiral g matrix 
reproduces the experimental data with no adjustable param¬ 
eter for the present scattering. Therefore, we can investigate 
chiral-3NF effects on proton and ^He scattering clearly. 

In Sec. II, we recapitulate the BHF method for the symmet¬ 
ric nuclear matter with 2NFH-3NF and the folding model for 
proton and ^He scattering. In Sec. Ill, the results of the fold¬ 
ing model with Gaussian chiral g matrix are shown for proton 
and "‘He scattering. Section IV is devoted to a summary. 


II. THEORETICAL ERAMEWORK 

A. Nuclear-matter calculations for 3NE 

We recapitulate the BHF method for the case of 2NFH-3NF 
[9]. The 3NF V 123 is hard to treat even in infinite matter. We 
then derive an effective 2NF Vi 2 ( 3 ) from V 123 by averaging it 
over the third nucleon in the Fermi sea. After this approxima¬ 
tion, the potential energy is reduced to 

i ^ {kik2\Vi2\kik2) A 

kik2 

+ ^ X! {^i^2k3\Vi23\kik2k3) A 
kik2k3 

= i^(fclfe 2 |yif|fclfc 2 )A ( 1 ) 

k\k2 

with the effective 2NF 

W? = ^12 + g^l2(3)) (2) 

where the symbol A means the antisymmetrization and ki de¬ 
notes quantum numbers of the z-th nucleon; note the factor 
1/3 in front of Vi 2 ( 3 ) in Eq. (2). The g matrix 512 is then 
obtained by solving the equation 

512 = V 12 + ^1*2*^0512 (3) 

for gi 2 with the nucleon propagator Gq including the Pauli 
exclusion operator. Here the single-particle energy for a 
nucleon with momentum k in the denominator of Gq is ob¬ 
tained by [9] 

ek = {k\T\k)+Re[U{k)] (4) 

with the single-particle potential 

fep ^ 

G{k) — \gx2 + gEi2(3)(l + Go(/i2)|fcfe )a^ (5) 

k' 
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FIG. 1: (Color online) Differential cross sections for neutron-proton 
scattering at Sin — 65 MeV in free space. The solid (dashed) curves 
represent the results of original (Gaussian) chiral t matrix. Experi¬ 
mental data are taken from Refs. [38, 39]. 


where T is a kinetic-energy operator of nucleon with the mass 
m and k is related to the incident energy E-m as Ein = 
{hkY/{2m) + Re[ZT]. When E{^ > 0, the single-particle po¬ 
tential is nothing but an optical potential of an extra nucleon in 
nuclear matter. Similar calculations, but in the second-order 
perturbation, of the optical potential in the framework of Ch- 
EFT was reported by Holt et al. [34]. The present formulation 
is consistent with theirs in virtue of the factor 1/6 in Eq. (5) 
[9]. 

In the present BHF calculation, the cutoff energy A = 
550 MeV is used with the form factor exp{—(q'/A)® — 
(q/A)®} both for V 12 and Vi 2 ( 3 ). The low-energy constants 
of chiral forces are taken from Ref. [35] as { 01 , 03 , 04 ) = 
(—0.81, —3.4, 3.4) in units of GeV“^, and the other constants 
{od,ce) = (—4.381,—1.126) are from Ref. [8]. Other sets 
of low-energy constants present in literature [36] are expected 
to give essentially same results. Furthermore, the variation of 
g matrices is much reduced in the effective 2NF level when 
3NFs are incorporated consistently [9]. In addition, the net 
effect of CD and ce is, small, when cd — 4:Ce- This rela¬ 
tion is well satisfied in various calculations for light nuclei 
in Ref. [37] and also for nuclear matter in Ref. [8] and the 
present work. As for U, our results are similar to those of 
second-order perturbation calculations [34] for the real part, 
but for the imaginary part the former is more absorptive than 
the latter. This may be originated in the full ladder-summation 
in p-matrix calculations. 

Figure 1 shows differential cross sections for neutron- 
proton scattering at E^n — 65 MeV in free space (in the limit 
of p — 0), where Ein stands for an incident energy in the 
laboratory system. The solid and dashed lines denote the re¬ 
sults of original and Gaussian chiral t matrices, respectively; 
note that the g matrix is reduced to the t matrix in the limit of 
p = 0. Thus the Gaussian t matrix well reproduces the result 
of original chiral t matrix. 

The g matrix can be classified with S, T, Ein and kp as 
g^"'"{kp, Ein)■ Hence U can be decomposed into lA = 
X;sr(2S' + 1)(2T -F 1)W^'^ with defined by Eq. (5) 



kp [fm'^] 

FIG. 2: (Color online) fcp dependence of at Em = 65 MeV for 
(a)^E, (b) ®E, (c) ^O, and (d) ®0. Squares (circles) mean the results 
of original chiral g matrix with (without) 3NE. The solid (dashed) 
lines represent the results of Gaussian chiral g matrix with (without) 
3NF. Eor ®0, the imaginary part is shifted down by 10 MeV. 

in which g and Vi 2 ( 3 ) are replaced by and respec¬ 

tively. Thus, 77'^^ means the single-particle potential in each 
spin-isospin channel. 

Figure 2 shows kp dependence of The squares and 
circles stand for the results of original chiral g matrix with 
and without chiral 3NF, respectively. The difference between 
the two results mainly stems from the 27r-exchange diagram 
in chiral 3NF. Particularly for nucleon and "^He scattering, the 
region kp < 1.35 fm“^ (p < po) is important. For the spin- 
triplet channels (^E and ^O), the 27r-exchange 3NF enhances 
tensor correlations and makes transitions between different 
states stronger, and eventually it makes the imaginary part of 
more absorptive. For ^E, chiral-3NF effects are large 
and repulsive, which corresponds to the suppression of A iso¬ 
bar excitations in nuclear medium in a conventional picture. 
The solid and dashed lines correspond to the results of Gaus¬ 
sian chiral g matrix with and without chiral 3NF, respectively. 
The Gaussian chiral g matrix well reproduces the results of 
the original chiral g matrix. The potentials in the parity-odd 
channels, and ^O, are small. 


B. Folding model 

The formulation of SF and DF models is summarized in 
Ref. [28], together with the relation between the two models. 
The folding potential is composed of the direct and knock-on 
exchange components. The latter component makes the po- 
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tential nonlocal, but it can be localized with high accuracy by 
the Brieva-Rook (local momentum) approximation [15]. The 
reliability of this approximation is shown in Refs. [40, 41]. 
The resultant folding potential U{R) is a function of the rel¬ 
ative coordinate R of projectile and target. The odd (^O and 
^O) channels of are almost canceled between the direct 
and knock-on exchange components, and hence U{R) is de¬ 
termined mainly by the even (^E and ^E) g matrices; e.g., see 
Refs. [24, 25]. The S matrices for NA and ^He elastic scatter¬ 
ing are obtained by solving a one-body Schrodinger equation 
with U{R). 

Eor ^°Ca, ®®Ni and ^°®Pb targets, the matter densities are 
calculated with the spherical Hartree-Eock (HE) method using 
the Gogny-DlS interaction [42]. The spurious center-of-mass 
(c.m.) motions are removed in a standard prescription [21]. 
Eor ^He, we take phenomenological density determined from 
electron scattering [43] in which the finite-size effect of proton 
charge is unfolded by using a standard procedure [44]. 


III. RESULTS 

Eirst, we consider proton elastic scattering at £^in = 65 
MeV from ^°Ca, ®®Ni and ^°®Pb targets. In Pig. 3, differen¬ 
tial cross sections da/dSl and vector analyzing powers Ay are 
plotted as a function of scattering angle 9c.m. in the c.m. sys¬ 
tem. The solid and dashed lines stand for the results of chiral 
g matrix with and without 3NP effects, respectively. Chiral- 
3NP effects are small at forward and middle angles where the 
experimental data [45] are available, since the scattering is 
governed by the potential in the surface region where 3NP ef¬ 
fects are small because of low density. Only an exception is 
Ay at flc.m. 60° for ^OCa and '^®Ni targets. Chiral-3NP ef¬ 
fects enhance the spin-orbit part of (7 (i?) by a factor of about 
30 %, which may be the reason for this improvement. We 
confirmed that chiral-3NP effects are small also for ctr. 

Next, we show the angular distribution of da/dil for "‘He 
scattering at 72MeV/nucleon from ®®Ni and ^°®Pb targets in 
Pig. 4. The solid and dashed lines denote the results of chiral 
g matrix with and without 3NP effects, respectively. Eor both 
targets, chiral-3NP effects are sizable at middle angles 0c.m. ^ 
20° where the experimental data [46] are available. 

The scattering amplitude can be decomposed into the near- 
and far-side components [47]. When a detector is set on the 
right-hand side of the target, the outgoing wave going through 
the right-hand (left-hand) side of the target is called the near¬ 
side (far-side) scattering. The near-side (far-side) compo¬ 
nent is mainly induced by repulsive Coulomb (attractive nu¬ 
clear) force, and in general the near-side (far-side) component 
dominates forward-angle (middle-angle) scattering. Eor both 
^®Ni and ^°®Pb targets, large oscillations seen in the range 
0c.m. = 5 - 20° are a consequence of the interference be¬ 
tween the near- and far-side components. When the scattering 
is dominated by the far-side component, da/dVl has no os¬ 
cillation and is sensitive to the change of nuclear force. The 
far-side dominance appears at 0c.m. > 20°. Chiral-3NP ef¬ 
fects thus appear in the far-side dominant angles sensitive to 
the change of nuclear force. 



FIG. 3: (Color online) Angular distribution of (a) differential cross 
sections and (b) vector analyzing powers for proton elastic scattering 
at 65 MeV. The solid (dashed) lines denote the results of chiral g 
matrix with (without) 3NF effects. Each cross section is multiplied 
by the factor shown in the figure, while each vector analyzing power 
is shifted up by the number shown in the figure. Experimental data 
are taken from Ref. [45]. 



FIG. 4: (Color online) Angular distribution of differential cross sec¬ 
tions for "‘He scattering at 72MeV/nucleon from ®®Ni and ^°®Pb 
targets. The solid (dashed) lines denote the results of chiral g ma¬ 
trix with (without) 3NF effects. Each cross section is multiplied by 
the factor shown in the figure. Experimental data are taken from 
Ref. [46]. 
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Figure 5 shows the central part Uce{R) of U for "^He + 
208pb scattering at 72MeV/nucleon. The solid and dashed 
lines correspond to the results of chiral g matrix with and with¬ 
out 3NF effects, respectively. Chiral 3NF, mainly in its the 27r- 
exchange diagram, makes Uge{R) less attractive and more 
absorptive. This repulsive effect of chiral 3NF in Uce{R) 
comes from the repulsion in the channel of The re¬ 
pulsive nature suppresses da/dft at 0c.m. > 20° for ^He scat¬ 
tering, whereas stronger absorption due to chiral 3NF better 
separates the far-side amplitude from the near-side one. 



R [fm] 

FIG. 5: (Color online) R dependence of the central part of the folding 
potential for "^He-l-^°®Pb elastic scattering at E = 72 MeV/nucleon. 
Panels (a) and (b) correspond to the real and imaginary parts of 
Uce{R), respectively. The solid (dashed) lines represent the results 
of chiral g matrix with (without) chiral 3NF. 


IV. SUMMARY 

We investigated the effects of chiral NNLO 3NF on proton 
scattering at 65 MeV and ^He scattering at 72 MeV/nucleon 


from heavier targets, using the standard BHF method and the 
^-matrix folding model. We evaluated the g matrix from 
N^LO 2NF plus NNLO 3NF for positive energy in nuclear 
matter. The same calculations for negative energies account 
well for the empirical saturation properties of symmetric nu¬ 
clear matter. Chiral-3NF effects are mainly originated in the 
27r-exchange diagram. The 3NF contribution in the chan¬ 
nel enhances tensor correlations to make the optical potential 
more absorptive. In the ^E channel, the 3NE effect yields a re¬ 
pulsion that may correspond to the Pauli suppression of isobar 
A excitation in the nuclear-matter medium in the conventional 
picture. 

Eollowing the Melbourne-groupprocedure [16, 32, 33], we 
provided the chiral g matrix with a 3-range Gaussian form by 
making a htting to the on-shell and near-on-shell parts of 
the original numerical g matrix. This Gaussian form makes 
the folding procedure much easier. The ^-matrix folding 
model with chiral 3NE reproduces the experimental data with 
no adjustable parameter for proton and ^He scattering. We 
found that chiral-3NE effects are small for proton scattering 
but sizable for "‘He scattering at middle angles 0c.m. ^ 20° 
where the experimental data are available. Chiral 3NE yields 
repulsive and absorptive corrections to UcEiR) for both pro¬ 
ton and ^He scattering. ‘‘He scattering is dominated by the far- 
side scattering amplitude at middle angles 0c.m. ^ 20°. The 
repulsive nature of chiral 3NE suppresses the far-side scatter¬ 
ing amplitude, whereas the absorptive nature of chiral 3NE 
better separates the far-side scattering from the near-side one. 
Chiral 3NE thus becomes sizable in the far-side dominant an¬ 
gle range. Note that the repulsive nature comes from the ‘E 
channel, whereas the absorptive nature from the ^E channel. 
Phenomenological 3NPs also make repulsive corrections to 
Uce{R) [17-19]. However, the origin of the repulsion is dif¬ 
ferent. It is interesting if the mechanism of producing the re¬ 
pulsive contributions is clarihed by analyzing scattering data. 
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